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An often ignored aspect of electromagnetic radiation from 
antennas is the characterization of their near-fields. A 
computer program, Numerical Electromagnetics Code (NEC), is 
validated for accurate near-field computations and applied to 
a model of a broadcast monopole. E- and H-fields are plotted 
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ieee KODUCT TON 


Peer 3 OBLEM ENVIRONMENT 
Maxwell's equations describe, for all space, the Electric 
(Band Magnetic (H) Fields that are the medium for trans- 


mitting information (virtually) instantaneously over 


extended distances. This practical application of electro- 
magnetic (EM) radiation has generated considerabie interest 
in its far-field* characteristics. The = and H vectors can 


be analytically determined when the current (I) on the 


generated) nas 


= 


Mama cor (trom which the field vectors are 


been described. For the far-field region of monopoles and 
dipoles, acceptable results are obtained by assuming a 
Sinusoidally varying current. However, when the field 
Weecors are Caiculated for field points at distances less 


miamera/2tT (the near-field region), significant errors occur 
between calculated values and empirically measured values. 
Mars ed2sicrepancy is mot the result of inaccuracies in 


~_ e-7 3 ? rey _ es c + — ie ° — a ~ Pa . 
Pemweli‘’S Equations, but can oe wholly attributed to an 


Pieecurasce desesipticn or the current. 


Mee farn-Trield, for moncrole or divole, is dcefined as cis- 
tances from the radi atior greater than ’4/27, wnere A 1s the 
wavelengtn. For free space, r=c/f where f=the Frequency at 
wmaeeh the nadiator is excited, and c=the speed of Light in 


— 
Le 
Tres space. 





in the past, the inaccuracies of the near field calcula- 
tions were not a problem--there was no practical application 
to which the energy could be applied. However, in recent 
years there has been an escalating interest in the close-in 
problem. This interest has been generated by an increasing 
proliferation of devices which use Radio Frequency (RF) energy 
in such a manner that the user is exposed to near-field 
radiation (e.g. microwave ovens, hand-held walky-talkies, 
Wireless telephones, RF sealers, etc.). 

The U.S. Navy was one of the first organizations to 
express an interest in these fields. The interest was the 
result of a concern for the hazards these fields might vose 
to fuels, ordnance, other EM operaticnal and test equipment, 
mieeecto personnel. The hazards posed by fuels, ordnance, 
Em@awotner EM equipment result in overt effects; the hazards 
posed to personnel are not so obvious. 

With virtually unanimous consensus, the literature 


Mdentiiies the heating of body tissue as the primary hazared 


W 


TO personnel. However, some distinctly non-thermal effect 
have been observed. Some of these effects [Ref. 14 include: 
a) minor changes in bicod properties 
b) a "buzz" heard by certain people when exposed to 
microwave radiation* 


ree Dues! 1s assumed to be 
Mere (ree) wather =5en the c 





c) abnormalities of the chromosome structure 

d) movement, orientation, and polarization of protein 

molecules 

e) changes in the transport rate at the blood-brain 

barrier 

f) comfort imbalances, such as epigastric distress, 

emotional upset, and nausea. 
The significance of these effects is not yet understood. 

The general public has not been exposed to near-field 

radiation hazards of sufficient intensity or duration to be 
cause for alarm; however, there are certain vocations in 
which the incumbents have an elevated probability of 


exposure to this type of field.* 


Peete lLS STATEMENT, SCOPE AND LIMITATIONS 
In order to provide more accurate data [Tor peopie 
researching biological effects of near-field radiation, this 


study will validate a computer program which computes the 


near electric, magnetic, and composite (peak) fields of an 
arbitrary radiator. The program will also be used to 
Compute these near fields for a model to be described beicw. 


The program, Numerical Electromagnetics Code (NEC), was 


develoned at Lawrence Livermore Laboratery under the 


“Sailors at sea or steevliejackse who effect antenna revairs, 
Hem instance. 





sponsorship of the Naval Ocean eViseemse Center (NOSC) and the 
Air Force Weapons Laboratory. It can be easily learned and 
applied directly by a researcher who has some knowledge of 
electromagnetic theory. 

The model to be used for the investigation will be a 75 
meter high broadcast monovole radiating above a perfectly 


Sem@awetine geround plane. The ground plane is located in the 


X-Y plane; tne monopole is co-axial with the Z-axis. The 


Hs 


antenna wili be excited at its base with a voltage sufficien 
to produce a radiated power of 1000 watts. The excitation is 
at a frequency of 1 Megaherts, which has a wavelength of 390 
meters; hence, the monopole is a quarter-wavelength radiator 
operating at resonance. A diagram of the model is shown in 
mroure 1. 

In simulating this model, NEC requires that the antenna 
be broken into virtual short straight segments. In 
@omsonance with this requirement, three different conftigura- 


tions will be used for computations; a model composed of 5, 


? 


Mmemend 25 sepments of 15, 5, and 3 meters, respectively. 
Using cylindrical coordinates, tne radial and parailel 

Peeetric fields, and the tangential magffetic Tield will be 

computed. Field values will be computed in .1l meter 


increments for points parallel to the Z-axis at distances of 


ci 


1 


ion oy 20, 100, and 2000 radii (1 radius .3 me S). 
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Field values will alsc be calculated along the radial axis 
Meeaxis) in increments of .5 meters for heights (h) of 2h, 
mm/s 2h, 1/4h, and 2/100ths h.* 

The near-field examined in this study is the cylindrical 
volume bounded by a cylinder of one wavelength radius 


centered on the Z-axis, rising from the ground plane to a 


height twice that of the monopole. 


* This is the approximate height of an average mans upper 
EOrsc . 
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weer PCN Or SELECTED ANALYSES 


feeeeCLASSICAL ANALYSIS 
The antenna utilized in this study is a member of the 


M_ptes1Centaitieda as 'thin-wire' antennas. The criteria 


ee 


which adumbrate the thin-wire approximation are [Ref. 2] 
a) Transverse currents are negligible relative to axial 
Sa preents 
b) Circumferential variations of the axial current are 

egligible 

c) The current can be represented by a filament current 
on the wire axis 
d) Boundary conditions on the electric field need be 
enforced only in the axial direction. 

By far, the largest subset of thin-wire antennas is a 
class known as short antennas. Short antennas are defined 
as having one dimension (h) much greater than the other 
dimensicns, but h is also much less than a wavelength, 
mypically n < i/8. Because of the predominance of this 
Sub-set, it is reasonable that it should be used as a model 
for classicai analysis. 

Based on Maxwells Equations and geometric considerations 


me shown in Fi 1 a. 


wv 


ure 2, the following non-Zero complex sc 


ey 
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field components can be derived: 
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where: n = the intrinsic impedance of the medium, I = antenna 


Mmerent, A = wavelength, 8 = 2n/rX = phase constant, h = 
monopole heignt, r = distance of field observation point 


From monopole base. Since there is no variation in $: 


The interpretation of these equations under classical 
analysis require two additional assumptions: 


~ 


aj} The distance to the field observation point is lLarce 
memariyesco the meonspole height, and 
bJ An harmonically time varying current is uniform over 


pe 


(D 


Beeant Or Ehe monopole. 
The first of these assumptions is made to focus attention to 
that volume of space where the lines of force have detached 


themselves from the antenna and are approaching the form ofr 


}-? 





a plane wave front. The second assumption is to forma 
First approximation of the current distribution on the wire 
Mmiteen 15 solvable in closed form. Extensive analysis over 
many years has shown this to be a valid assumption when the 
antenna is excited at a single point. 

In this study, we will be interested in fields where the 
observation points are not at distances considered large 
relative to the monopole height. When this is the case, the 
diagram of Figure 2 must be modified as shown in Figure 3%. 
From the geometry of this figure, the following non-zero 


complex scalar field components can be derived: 
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H. = (2 + e - 2{cosBh)e 
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*Fquations (1) through (3) are expressed in spherical coor- 
dinates while Equaticns (5) through (7) have been transformed 
me Cylindrical coordinates. 
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Percurprent of the form: 


mez) = T,singh mooie < 2 = 1 (3) 
has been inserted into Equations (5) through (7). 

Equations (5) through (7) indicate that each field 
component is the phasor sum of contributions from three point 
sources. The first source is located at the top of the mono- 
pole, a second image source is located in a symmetric position 
when the first source 1s reflected through the ground plane. 
The third source is located at the base of the monopole; its 
relative amplitude is modified by the (-2cosg@h) factor. 

Since we have a solution for computing the field points, 
it is now reasonable to determine the range over which it 1s 
valid. By inspection we note that for vr = \/2n in Equations 
Mm ehrough (3), the magnitude of each non-Zero term in the 
parenthesis is equal. This implies that r = \/2n may be a 
fmeem@ettion point. Further investigation indicates that for: 


S , : 2 
Bees \/ 21: che l/r terms in Ee and H,. and the l/r 


( 


eo bet 


term in E_ are dominant 
r 
3 bs = sam] 
mee C<\/27: the l/r terms in EE. and En and the l/r 
term in Hi. are dominant. 
Since we know that power drops off as the inverse square of 
distance, and that the power flow at a field point can be 


SPerined in terms of its E-field as: 


{) 


| 
or 





P=E S/n (9) 
2 2 
we can extrapolate from a) above (E 8/ne« k/r”. where k = a 
constant) that the far-field begins at r = A/2n and continues 


out to infinity, while the near-field begins at the monopole 
M@eameontinues out to r = A/2r.* 

It is interesting to note the complementary relationship 
that is maintained by the power in this transition regicn. 
Instantaneous power flow at a point in space is given by the 


Semolex Poynting vector: 


st: 
HW 


Wee eat) (10) 


For the short monopole, the complex Poynting vector is 


oriented in the radial direction and given by: 


NO 


-j8r)2/8 
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The real part of this equation gives the éverage rate of flow 
‘ ae <5 ; a. ae 
of Real power («l/r ) while the Imaginary part is proportional 


to the difference between average magnetic energy density 


Dlnis a>cument is Offsred not as a proof, but as a reascnable 


Pationale which can be empiricaily verified. 





: : fe) : 
and average electric energy density («l/r ). Interpretation 
Seeents equation indicates that for: 


: : iE 
a) r>i/2n: Reactive power flow decreases rapidly («-,), 


r 
Real power flow decreases not so rapidly (=>) 

r 
b) r = iA/2n: Average Real power flow and average Reactive 


power flow are equal 

ec) r<A/27: Reactive power flow becomes much greater 
5 : 2 

(<r) than Real power flow (er’). 


The transverse component ef the complex Poynting vector: 


)-cosesine(-j(8, + 5)) 6 (12) 
Q spe 


is completely Imaginary, indicating no Real power flow, and 
it falls off rapidly with increasing distance due to l/r? and 
Ue gem EPibutions.* 

The deviant behavior of the fields encountered at 
distances r<A/2n7 is the result of moving from a region where 
the lines of force ere closed on themselves and approaching a 
plane wave front, to a region where the lines of force still 
maintain a strongly spherical wave front and some of the lines 


are still attached to the antenna. Even in the simplest of 


“Also noted is the fact that the Electric field changes from 
Mmaptical polarizaticn to Linear polarization in passing 
me@emer>A/en tO value of v<ik/27 when 6 # 0°, 90°, or 180°. 
For these angles the pelarization of the electric field is 


[Meeear on bOotn Sides of r = A/27. 
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Gases, the current distribution in conductors is complex. 
This complexity 1s the result of non-uniform currents and 
charge densities, capacitance and inductance generated by 
these non-uniformities, and discontinuities from impressed 
moecineg functions. 

It is well established that as the field observation point 
moves within the \4/2n distance from the antenna, that the 
measured values deviate from those predicted by Equations (5) 
through (7). It is also well established that Maxwells 
eduations are valid for all space. Consequently, the con- 
elusion is that the sinusoidal current distribution is not 
Suerici1entiy accurate for determining field points within 
the near-field region. The problem then is to find a better 


s the domain 


} 


Meeeription of the current distribution. This 
of contemporary analysis. 

Before proceeding, it 1s necessary to reflect that the 
analysis performed thus far has been for a short antenna. 
The model for this study is a resonant monopole; when the 
antenna 1s lengthened to resonance some interesting 
characteristics appear. Primary among these 1s the fact 
that the reactive compcnent of the impedance vanishes, which 
says the input reactance goes to zero. The current on the 
antenna and the voltage required to place it there are 
determined solely by the resistance, the real part of the 
impedance. When a sinusoidal current is assumed on the 


antenna, the tangential field component is given by: 


20 
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The (-2cosgh) term representing radiation from the base of 
the monopole has disappeared. The monopole appears as a 
point source radiating from its height, and an image source 
found at the reflection through the ground plane. 

Also, the resonant condition occurs at a lengtn less 
than that expected from the A/4 relationship. This is the 
result of an "effective" lengthening of the antenna by 
Capacitance and "fringing" at the end-cap of the wire. 
Capacitance is caused by a build-up of charge at the end of 
Petre O- finite width. The change in direction of the 
geometry at the end of the wire also changes the directicn 
of the lines of force emanating perpendicular from it; this 
is the fringing effect. These factors cause resonance to 
occur for an antenna length of h = A/4 ~- 6 where § is a 
function of the capacitance caused by the two effects 
discussed abcve. 

The short monopole development discussed above is very 
restrictive and cannot be applied directly to resonant 
monopoles. However, when the resonant peculiar characteris- 
tics mentioned are taken into consideration, the resulting 
development will follow the same procedure as used for the 
short monopole. For the resonant menepoie development and 


-~ —_ cao! 
Mesulizs, see Ref. 4. 
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Pree CONTEMPORARY ANALYSIS 
Schelkunoff [Ref. 3: pp. 370-374] identifies three methods 
for obtaining a closer approximation to the actual current 
distribution on an antenna. They are: 
a) The Integral Equation method 
b) The Sweep-Off method 
c) The Mode Theory of Antennas method. 
Each are discussed below, in succession. 
1. The Integral Equation Method 
The assumption is made that the current is distributed 
Saetne suriace of a hollow cylinder which is divided into 
filaments of angular density (1(¢z')/2n) and angular width 
(d@'). The Kernel Function is modified to allow for 
variations in % and is expressed as an integral around the 


Memeeimierence of the wire as: 
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straight current filament: 








The field intensity for a 
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is modified by inclusion of Equation (14), with the boundary 


Semaation: eee ED = 0) enforced to give: 
i l "2320 (16) 
mE-(2) = f (Sey +8°%) T(zt) az! 
: SE ae 32 


This Integral Equation (from which the name cof the method 


derives) is the limit (as N+) of a series of equations of 


the form: 
2y 
: N SRN Gir 2) aren Ae) 
mez) = 5 =~ [ zi + 4 CZ) i ee Z >) 
Z Bee we a n 
roles te 0 Z 
Cl 
free 1s also of the form: 
aie 
wee(z) = Z se a Cie} 
Z, mn 


Equation (16) is a ‘circuit equation! for an antenna and can 
De solved, in a manner analogous to the use of Kirchhoffs 
Equations for lumped networks, to an arbitrary degree of 


Peeuracy through the use of Equation (17). 





Peis wiisurlecrtive at this point to pause and 
consider how accuratelyotrinitial approximation of the current 
Spewodehed reality. Integrating Equation (15) by parts and 
Beweelng ~E*(z)=0, as it 1s on the surface of the antenna, we 


eoeain: 


© iia 2 
(T(z")%(z,2") + Te2S) = op 2 (LE + 32108 ¢z,2") aztec 
az! 
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CaS). 


This equation shows that as the antenna radius approaches 
zero, the first term of Equation (19) approaches a constant 


iegeme for all z2, except 2 = Z, or 2 = Z,. Sur iseeently fer 


m@emethne ends this constant is small. The ‘(z,z') factor in 


maeesecond térm is infinite at z = z' and large in its 


2 

= 
—7 ar 3*T) factor must be small. 
az 


Therefore, as the radius of the antenna approaches zero, 





Meeinity; hence, the (¢ 


om Z ta eo yan 

576 I) approaches zero, which implies I(z) approaches 
az! 
a Sinusoidal form. 


2. The Sweep-Off Method 





This method starts with Equation (15), the tangential 
field equation for a filament current, being integrated twice 
Byeparts. Then an approximation to the current distribution 
on a transmitting antenna of finite radius, assumed to be 


Seven ly: 


24 





Mey =f Cosez' 9? B sin8gz! (20) 


is substituted into the integrated form of Equation (15) to 


yield: 
. =F, =387, eu. 
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This equation is an exact expression for the electric field 
intensity parallel to a sinusoidal current filament. The 


antenna is then subjected to a compensating impressed field 


Y 
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ae : 2 : : 
Semana une correction terms from the first approximation 
are computed. The name for this method resuits from the 


concept of 'sweeping-off' the surface of the antenna the 
residual tangential electric intensity by applying an equal 
and opposite intensity. 

It should be noted here that although Equation (21) is an 
exact expression for a thaw ag does not vanish along the 


Ty 


Filament. This apparent dichotomy results from the erroneous 


Sewer On Of UNITOrm Convergence for E_. In actuality, the 
convergence of E. 1s non-uniform at o = 0, where o is the 
distance from the filament. For an antenna radius r = a, 





there exists a remainder term T, (2) nGipcte Ginseeom wes 1 (Zz): 
I, (2) Bomsutiacienteto cause E,(a) = ees sc Ole T,(2)+0; 
Gemsequently, its contribution to E. for any fixed value of 
9, greater than a, will approach zero with a. In the limit, 
E. memowven by Equation (21) for 5>0, but at p = 0 Equation 
Seeedoes not hold. 
3. Mode Theory of Antennas 

This method is based on the solution orf Maxwells 
Equations subject to boundary conditions at the antenna 
surface, and the surface of the source of power. The method 
Memeasts Or Calculating modes of propagation consistent with 
the boundary conditions at the lateral surfaces of tne 
antenna, then the modes consistent with free-space prcpaga- 
tion. These modes are combined to satisfy all boundary 
conditions and the current is obtained as the sum of two 
components: the TEM cr Principal wave, and the Complementary 
wave ccnsisting of all higher order modes. The principal 
mode is governed by ordinary transmission line equations 
Men cerms of distributed capacitance and inductance). 

+. Discussion 

Of the metheds discussed abcve, 'Sweeping-Off' is the 

most elementary and consequently, the least precise. The 
Mode Theory of Antennas excites a current in the conductor 
for the TEM mode, and a current for the complementary wave, 


which is a sum of currents, one for eacn of the higher order 


modes excited. The excited modes are strongiy dependent cn 





the geometrical configuration of the system; hence, a 
generalized solution requires a number of terms which rapidly 
approaches infinity. This overwhelming problem of accounting 
for all terms limits the attractiveness of this method. 

The Integral Equation method apnears most promising 
for generating a better approximation to the current. Many 
different approaches to solving the equation have been 
proposed and virtually all of them lend themselves to 
sOluticn by digital computer. Three of these approaches will 


meme rscussed in the next section. 


C. NUMERICAL ANALYSTS 
1. Hallen's Approach 
One cf the earliest solutions to the Intesral Equation 

was the approach proposed by Hallen,* the basis of which is 
successive approximations to the current. Hallen uses a 
Meqeit ional Green's Function along with a tacit assumption 

of a homogeneously distributed current** over the conductor 
pmmrace. the 'E' Tield inside the conductor (expressed in 


terms of current and skin effect resistance) is equated to the 


“Hallen's approach is relatively long; the details will not 
be presented here. For an outline of the method and a discus- 
Sion of the important steps and results, see [Ref. 5]. 


“tHallen's tacit assumption is that the impressed voltage is 
distributed over a finite length [Ref. 3] (which, for a homo- 
geneous conductor, implies a homogeneously distributed current). 
/Meetne stated assumption of a driving cotential discontinuity 

Mee toa, eugorced, Hallen'’s method breaks down [Ibid.]. 
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fE* field outside the conductor (expressed in terms of 
vector potential, 'A'). The wave equation is solved as a sum 


Gea complimentary function and a particular integral. A 


g 
@enstant, Co» in the solution is evaluated in terms of input 
conditions at the terminals, and the vector potential, A, is 
expressed in terms of the antenna current. C. and A are 
inserted in the solution for the wave equation obtaining an 
Mmeegral Equation in current, which is Hallen's Equation. A 
partial solution for the current is obtained by evaluating 
one of the integrals so the current is expressed as the sum 
of several terms, some of which also include the current, I. 
Neglecting certain terms in I, an approximate (zero order) 
solution is obtained for I. This value is substituted back 
in the current equation obtaining a first-order approximation. 
This process is continued to yield successively higher orders 
of solution. The process is stopped when a sufficient degree 
@eeacecuracy is obtained, the constant Cy (the coefficient for 
the quadrature term of current) is evaluated and an 
Peempeocic expansion for the current is obtained. 
2. Parabolic Interpvolative Fit Approach 
A more recent approach (1975) is one proposed by 
Chang et al, [!Ref. 6]. Chang begins by assuming a4 surface 
Berrent on the antenna (as opposed to a line current 
approximation) and defines a new Kernel (Green's) Function 
to account for axial variations in current. He also derives 
< 


peolsserent expression for E_ from the relationship: 
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B= -j(0(0+R) + 3°R) (22) 


fem@en he Claims 1S more accurate, since centinuity of higher 
order derivatives of 'I' do not have to be specifically 
considered at segment boundaries when solving the Integral 
Equation. He then assumes a Parabolic Interpolative Fit for 
the current distribution over each segment (the distribution 
over the whole antenna is not necessarily parabolic) where 
precise values of current are exactly defined at sample points. 
He also defines three moment functions that describe the 
Kernel Function and its variation with respect to # and 2 

over a segment. The non-zero field points are then expressed 
as double summations of the moment functions over each inter- 
val. To avoid the double summations, he splits the Kernel 
Bem@erton into two factors (a technique used by Hallen, as 
weli) which facilitates the computation of the dominant term 
(the result of which is a complete elliptic integral of the 
first kind), then redefines the Kernel as an elliptic integral 
plus an approximate remainder. With this Kernel he then 
computes the moment functions defined above wnile redefining 

a remainder term for three different special situations. 

These inoment functions can now be evaluated numerically 
through adaptive integration techniques, such as Gauss-Legendre 
Quadrature, by splitting the integrand functions into Real and 


Imaginary parts. 
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3. Bessel. Function Fit Approach 

Another recent approach (1981) is that proposed by 
Balzano et al, [Ref. 7]. Balzano also begins by redefining 
m@eekerne! Function to account for azimuthal variations, but 
Sees tt into the form of an integral of a Hankel Function of 
the second kind, which also can be represented as an infinite 
Summation of Bessel Functions. The magnetic vector potential 
is then expressed as a double integration of the current 
Semsity on the conductor times the integral or the infinite 
Summation of Bessel Functions described above. In this form, 
all that is needed is an accurate description of the axial 
Current distribution, I(z') over the cylindrical antenna. 

From the relationship of the vector potential with 
the E and H fields, an expression can be written for E_ as a 


memetion of the integral of I(z): 
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I¢(z') is expressed by a suitable linear combination of 
orthonormal functions which satisfy tne boundary conditions 
for the current. Then the current expressed as an infinite 


Miwon cnese Ortionormal functions is equated to the Fourier 
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Heasrorm Of the current expressed as a function of the 
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Wector potential: 
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where: SNS! 1s the Fourier Transform of the current 
harmonics, oe is an infinite series of constants. Equation 
Mmeeyeis substituted into Equation (23) and solved by applying 
Smerecins Method.* This yields an infinite system of 
equations which can be truncated and solved for the unknown 
Ave A similar treatment can be apolied to the other non-zero 
field equations to cast them into solvable form. However, 
the equations so defined converge slowly; therefore, addi- 
tional manipdulation is required. The additional manipulation 
becomes rather involved and will not be discussed here. 
Mee@ails can be found in Ref. 7. 

4Y. Discussion 

All three of the above methods are theoretically 


Teakovghe 
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sound, all three give reasonable results cover large bD 
of the near~field region, and all three lend themselves to 
fmerteal solution by computer. Of concern here is the 
amount of valuable computer resources (time, core, and cost) 


° 


consumed in making these calculations. Balzano indicates the 


*Galerkin's Method solves an Integral Equation via the Method 
6: Moments to be discussed below. Specifically, Galerkin's 
Method sets the ‘leighting Functions equal to the 3asis 

React Ons. 
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expressions he derives for field components are slow in 
converging, and even after modifications to his basic 
derivation to hasten the convergence, he states that more 

than 1500 harmonics (a 1500th order linear system of equations) 
are needed to satisfactorily match boundary conditions. 

Chang does not indicate directly the amount of 
computer resources consumed by his method, but he does imply 
his method is adopted because a simpler method (The Method of 
Moments) 1s inaccurate: 

eee Chie mumerical computation, the integral from ~-h to h 
is usually subdivided into a finite number of segments; in 
each segment the current 1S interpoiated by a polynomial 
with coefficients expressible in terms of values of current 
at a few sample points within that segment. The higher 
derivatives of I(z) are therefore discontinuous at boundary 
points between any two segments." [Ref. 6] 

One of the intentions of this study is to show that 
the Numerical Electromagnetics Code (NEC), which uses the 
Method of Moments, has been modified to account for the 
discontinuities of the first derivative at segment boundaries, 


and that this modification is sufficient to yield results that 


are accurate for engineering analysis. 
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Wee, Numerical Electromagnetics Code, 1s a computer pro- 


gram for the analysis of the electromagnetic response of 
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Ca 
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emeenhnas and other metal structures. It computes a numer 
solution to integral equations that describe the currents 
induced on a structure by voltage or current generators and/or 
incident fields. It is the latest in a series of programs. 
The rirst, called BRACT, was deveioped by MBAssociates under 
funding from the Air Force Space and Missile Svetems Command. 
This was followed by AMP, again developed by MBAssociates, 
this time under an Office of Naval Research contract. NEC 
evolved from AMP and was developed by Lawrence Livermore 
Laboratories with funding from the Naval Electronics Systems 


Command. It is a user friendly ccde which incorperates the 


moullowing: 
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paeeemecoma!l EQuation €1,£.) for modeling the current on 
smooth surfaces is combined with an I.E. for surrent modeling 


Oman wires to describe the electromagnetic response of an 


C2) 


arbitrary structure. The structure, or parts of it, may be 


active or passive, located over a perfect or imperfect ground 
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plane, end may have lumped-element or distributed loading. 
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xCLtation may be 
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structure, an incident plane wave with either linear or 
elliptic polarization, or the field may be due to a Hertzian 
dipole. Output may include current and charge density, power 
gain or directive gain, near- or far-zone electric or magnetic 
fields, impedance or admittance, total radiated power or input 
power. it is suitable for either antenna analysis or cross- 
section scattering and electromagnetic pulse (EMP) studies. 
The code wili handle junctions of wires of uneven radius, free 
ends of wires that have finite radius, wires of variable radius, 
and coupling between wires. 

NEC utilizes the Gauss-Doolittle Method for solving the 
matrix equations generated by the Method of Moments when 
Solving the integral equations. It also allows for use of 
rotational or plane symmetry to reduce computation time. When 
the impedance matrix is too large to be contained in core. NEC 
has the option of storing portions out of core. It ailiows the 
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Mself—interaction matrix' for a structure to be computed, 
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factored for solution, and stored o or 
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Porm aA new antenna that enters this environment requires only 
mmeweyalliataon of the 'seli-interaction matrix' for the 
antenna, the mutual antenna to environment interactions, and 


Metrix Manipulations for a partitional matrix soiution. 


Pree ZONING CONSIDERATIONS 
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NEC is a discrete sampling code where 
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Bmates) ©O which the Electric Field Integral Equation (EFIE) 
or Magnetic Field Integral Equation (MFIE) are applied. As 
with any description of the real world, there are aporoxima- 
mens, buc because of the versatility in modeling the geo- 
metry of a structure, its avproximations more closely 
resemble nature. The resemblance is strongly influenced by 
the choice of zoning (i.e. dissecting) the structure in the 
MmeegPam. Ihe smalier the geometric elements, the closer the 
model comes to reality. However, the smaller the elements, 
the larger the number of elements, which means the larger the 
Matrix Of equations and hence, the more costly the solution. 
There is a point beyond which smaller zones do not yield a 
substantially more accurate solution; it may be necessary to 
rerun the problem with increasingly smaller elements to find 
the point of diminishing returns. The choice of vroper zoning 
then is gained with experience and becomes as much of an art 
as it 183 a science. The guidelines for the science asvect 


ee as foliows. 


Segments should follow the paths of conductors using 
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Diece-wise linear fit on curves. Generally, segment lengths 


i) 
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(A) should be less than .14; shorter segments (.05\) or ie 
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ag 


may be needed at critical regions (junctions or curves 
: -2 . ae 
segments smaller than 10 A should be avcided since the 


Samiiaricy or constant and cosine components lead to 


~” 


numerical inaccuracy. The radius of the wire (a) relative to 
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me@epends on the Kernel used in the Integral Equation. Two 
options exist. The thin-wire Kernel models a filament current, 
while the extended Kernel models 4& uniform current distribu- 
tion around the segment surface. The field of the distribu- 
ted current 1s approximated by the first two terms in a 

° ° ° 2 2 
Series expansion of the exact field, in powers of a’. The 
my I ee ee» oe 
first (a ) term is identical to the thin-wire Kernel; the 
Peeone t2rm extends the accuracy ror larger values of a. 
Both Kernels incorporate the thin-wire approximations (see 
Seerton 1f.A) and both require 2ra/A<<1. The thin-wire 
Kernel requires a A/a>3; the extended Kernel relaxes this 
to A/a>2, These values ensure errors are less than 1%. The 


extended Kernel is used at free wire ends and between parallel, 


connected segments. The thin-wire Kernel is always used at 
bends. 
poe SUrtaeces 


A eonducting surface is modeled by small flat surface 


i-hy 


Patches which conform as closely as possible to curved surfaces. 
The parameter defining a vatch is a normal unit vector, 
criginating from the center of the patch, defined in Cartesian 
coordinates. Each patch must be connected by a wire at the 
Hatch center for the program to integrate the surface current. 
the code divides each patch into four equal patches about the 
feeeemeta, along the unic vector lines describing the surface 
Peete Sa-eh. An inteYrpolaticn Tunction is apo 


rour patches to r 
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Meeenes, and the function is numerically integrated. VYatches 
with wires connected to them (active patches) should be 
chosen to be approximately square with sides parallel to the 
unit vectors defining the surface. Only one wire may connect 
to a patch, that wire may not be connected to another patch, 
and it may not lie in the plane of the patch. A minimum of 
about 25 patches should be used per wavelength of surface 
erca;, the maximum size for an individual patch is about .04 
Square wavelengths. The number of patches used increases, 
and the size of each patch decreases, as the radius of 
Sm@evature decreases. Smaller patches should be used at 

edges since the current amplitude may varv rapidly in this 
region. Long narrow patches should be avoided. Patches are 
restricted to modeling voluminous bodies with closed surfaces; 
parailel surfaces on opposite sides cannot be too close 
fecether. 


me 6 6GGPOuUnGd «~ Plare 
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or a perfectly conducting ground, the code generates 
meee tected image. Structures may be close to, or contact, 
mee eround; however, for a horizontal wire: 
where a = wire radius, h = height of wire axis above the 
ground plane, and h/a>3. 


A finitely conducting ground may be modeled by an 
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image modified by the Fresnel Plane-wave Re 


Meemrte©ieies. il0h2zs method is fast, but of limited accuracy 
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Pee sneuld NOT De used For structures close to the ground, or 
— 2 


08) 





having a large horizontal extent over the ground. The 
Sommerfeld/Norton model uses the exact solution and is 
accurate close to the ground; the horizontal restricticn is 
the same as for a perfect ground. This model is only used 

for wire-wire interactions, for surfaces the code reverts to 
Fresnel Reflection coefficients. At the present time, a 
ground stake cannot be used, but wires may end on a perfectly 
conducting ground plane if the derivative of the current at 
the end of the wire is zero. For wires, there are options for 
a radial-wire ground screen approximation and a two-medium 


ground approximation based on modified reflection coefficients. 
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These options allow considerable savings in computational time 


when limited accuracy can be tolerated. 


See eNTEGRAL EQUATIONS 
eeetectric Field Integral Equation (CEFIE) 
The EFIE is used for thin-wire structures with small 
@onductor volume. It is derived from the electric field 
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representation for a current distribution confined to 


Surface of a perfectiy conducting body: 





Miers: £' is the source point, Y is the observation point, 
feeeets 2 Surface current densizy, Glr,r') = (8 IT + vy)elr.r'), 
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impedance of medium. The principal value integral #. is 
indicated so mr’ in the limit. 
An integral equation for the current induced on S by 
2a 


an incident field E” can be obtained from Equation (25) and 


the boundary condition: 


Mee iEawc) + E> (s)] = 0 Oe 
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where reS, n(r) is the unit normal vector at Yr, 1s the 
Field due to the induced current ee This integral equation 


1s given by: 
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This vector integral equation is reduced to a scalar integral 
equation for the thin-wire model. With the surface current 
density replaced by a filament current and the boundary 


condition enforced in the axial direction, this gives: 
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BecrG) = Soy nee Gee een enor )ds' (28) 
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mgeme: S = distance along the wire axis at r, and s' = unit 
Mm eerom along the wire axis. This is the EFIE solved by NEC. 


ee eo e facid tnceseal Equation (MFIE} 


The MFIE is derived from the integral representation 


meretne Magnetic field of a surface current distribution J. 
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et J. is induced by an external incident field 3”, chen the 


total magnetic field inside the perfectly conducting surface 


must be: 
A*(r) + AScr) = 0 (30) 


si e e ° ° —S ° 7 
mowee- 1 as the incident field, and H is the scattered 
meld. Let rr, along n(r_) from inside the surface, then 


Equation (30) with Equation (29) substituted becomes: 
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The limit is evaluated using a result of potential theory to 


yield: 


LO 








MEG eal oO in) ee nth) 
'@) © S cS TT S © 
mice’) xX —— e(e on") 1] dA’ (32) 
S or O 


This vector is resolved into two scalar equations along the 


orthogonal surface vectors Th and Lr where: 


t,(p_)-H4(p_) = -1/2 t,(P_)-T_(e_) 
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NEC models thin wires using EFIE and for surfaces 


uses MFIE. 


For a structure with both wires and surfaces, nr 


in Equation (28) is restricted to the wires, with the integral 


for E°(r) extending over the complete structure. 


Bomation (238) 
Meiacion (27) 
restricted to 
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wires, while the more general form, 
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Pauamiee soy and (35), wlth the 
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H- (r) extending over the complete structure. 


coupled integral equations are: 








surfaces: 
N Pa a 
ae 7 {\ 4 seh 
| r= ra i Ce (ese = a 
is 
; n on ees 
anes: f Ee, [8 oe 


S 


ut 


excluding wires: 


5 


= Saya aoe ds 


) 


! 


Sas 


cm r.r') dA! 
Sa BS ok ) dA 


= = t,(p)- ¢ I(s')(st x a ooo a ds * 
2 7 je 
1 Cn Cas i; es 
- 1/2 ct, (xr) Co Gay, =e ee io ee. 
abe = rar S Ze 
al 
- (FCS) X se, geet) ] da er 





and 


I(s')€s! x =r g(r,e')] ds! 


14) 


Ty 
fii. 


an am fe A as : 
2 Cae re t, (9) ‘ 


“1/2t,(B) +s Ce) 6 Be t, (B) 


t 
o- 


U 
al 
4 
ee 

te 


- (5_(e") X gop g(@,7")] dal (38) 


Sees HiOD OF MOMENTS AND THE MATCH POINT 


The Method of Moments is a techniaue whereby an integral 


equation is reduced to a system of linear algebraic equations 
which are easily handled by high speed digital computers. 


1. Mathematical Concept 


The method applies to an inhomogeneous linear 


weeracor of the form: 
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linear operator with domain D, containing the 
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function t, which is an unknown response to an excitaticn e, 


which is given, and found in the range cf L. 


UQ 


BeeomoOlicton = 2O sGUarion (39) exists and is unique 


meeeall! @€, then the inverse operator, L ~~, exists such that: 
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Peet 2Seoresents a solution to Equation (39). 
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The unknown function, f, may be expanded in a series of 


N linearly independent Basis Functions spanning D 


Peete ituling Equation (41) into 


the inner product with a set of 


weighting functions {wi}, defined 
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linearly 
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eo. and taking 


independent 


range L and spanning 


sub-space 3 results in a set of equations for the 


coefficients, Age Om povation (4d). 
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momeeDe written in matrix notation as: 
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This set of equations: 
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Pet ee et aoe Nn Monee esaUatzen (43) exists and is 
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Peace top all e, then the inverse operator. L ~, also 
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peemem 25 a solution to Equation (43). 

The efficiency of computations and accuracy cf 
solution is largely dependent on the choice of Basis Function. 
Factors which should guide this choice are: 

a) accuracy of solution desired 

b) ease of evaluation of matrix elements 

c) size of the matrix that can be inverted 

d) the realization of a 'well-conditioned' matrix. 
There are two general classes of Basis Functions from which 
to choose, entire domain and sub-domain. The sub-domain class 
has fewer elements; hence, its execution time is less, it 
Simplifies the evaluation of the inner product integral and 
still ensures the matrix [G] will be 'well-conditioned'. 

PeeecC's Application 

Nees tinterpoletion Of Current On wires enforces 
continuity of current and linear charge density at segment 
ends. This imposes additional restrictions on the current 
and 1ts derivative with respect to distance along the 
memaweror, NEC forces the current on tne conductor to match 
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che integral equation at specific points; these 'match po 
are located at the center of each segment. 


Om Of the Method of Moments starts 
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NEC's applicat 
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Peeemechne Inteeral Ecuations: eEquation's (23) and (34) te 
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me), Or Equation's (36) to (38); the equation selected 
depending on the structure being modeled. These integral 
equations are the basis from which the system of linear 
equations will be generated. 

A source model with a known excitation voltage can 
be related to the E-field impressed on the antenna; hence, 
a (= e from Equation (39)) is known. We know we are locking 


Berea more accurate description of the current distribution 


mone diatom (39)? 
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@meche antenna; hence, the current Tf (= f 
1s the unknown response. The operator L (an integral 
operator in this case) is the composite of all other terms/ 
Factors (other than E and I) found in the integral equations. 


Hence, these terms can be cast into the form of Equation (39) 


To generate the system linear equations, we have 
divided the structure to be modeled into a number of segments 


according to the zoning guidelines discussed above. The 
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current distribution on the antenna is now approximated a 


the summation of currents found cn each segment. That is: 
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which has the same form as Equation (41). Here g. are the 


Oo 


Beeis Functions and Le is the system of coefficients. For 


ares : 


iieea= ye + B. Sin€g(s=-s.) + C. cos(s(ses.))* (47) 
1a i) 1 J J 1 
where: 2s = center of the jth sub-section, and 8 = free 
eSpace wave number. As a result of current continuity 


considerations at segment ends, each sub-section will have 


Pity one unknown. 


For patches, the surface current density ls: 


Meer) =u kK. +v K . C43) 
S We] Vi 
where: Ae and Ro are vector components or the surface 
current density, Wee); on each patch. Each K. is constant 


Over the patch, which is usually rectangular, but can be any 

reasonable shape. Again, current continuity considerations 

for the patch yield only two unknowns, one for each direction. 
Because of the nature of the Kernel functions in the 


integral equations, the choice of basis functions is more 


Sete 1cal on wires than on surfaces. 


“This form has been found to provide rapid 
gence and fields of sinusoidal currents are 
zn closed form. 


solution conver- 
easily evaluated 





Zarolplidedigey. the Testeof thisusecticn will be 
Meetrictred to the case for only wires, as this is in 
Semoonance with the system for this study. Assuming the 
Wire 1s aligned with the Z-axis (as it is in this study), 


Equation (45) can be re-written as: 


i 
Pez = [, £ (zt) (49) 
Mimere the z' indicates the current is on the surface of the 


wire. Equation (46) can be rewritten as: 


Mem = 7 iL. g.(z") CaO) 
= ce! 


Fquation (50) defines the set of sub-domain basis functions 
memoe Used in computations. For finite N, the sum of fs Pez) 
cannot exactly equal the general current distribution*, so 

the = ae need to be chosen as close to the actual 
Gistribution as possible. This is accomplished by restricting 
maessupport tor Le gue? to a localized regicn of the surface 
near the center of the segment. Fach basis function has a 


peak on the segment on which it 1s defined and extends only 


*The approximation can be brought arbitrarily close to the 
Peeqetesolucion by choosing N sufficiently large. Fortunately, 
reasonable accuracy can be obtained, in most cases, witha 
mematively small number of segments (N 2 10/)). 


tes 
CO 





onto segments to which it is directly connected, going to 
zero with zero derivative at the outer ends of connected 
segments. 

Equation (50) is then substituted into Equation (49) 
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where: I, and the summation have been moved outside the 
J 

domain of L as a result of its linearity. 


The next task is to identify a set of weighting 
| 


monet lOns , twit. When the weighting functions are set eaual 
1 
Galerkins Methed.* In NEC, We is chosen as a set of Dirac 


| to the basis functions (l.e. w. = Bs); this is known as 
emecra Functions: 
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Bae-ee- = is the set of match points found at the center of 
eacn segment Ai. Use of this weighting function resuits in 
a point sampling of the integral equations known as the 


Peoellocation method" of solution. 
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ae ae ae ee 


Memploved =n Hallen's approach (See Section II.C.1) 








Vehngeen is swerehting function (Equation (52)) and 


imieimes the inner product of each side of Equation (51) gives: 


N 
eee ee OL. <w., Lie (z')]> Coes) 2 
aL he 5 iL 
jel 
which has the form: 
my = (1) (72) (54) 
where: [Vl] is a matrix of weighting functicnsapplied to the 


Known impressed field, [IJ is a matrix of unknown coefficients, 
and {Z] is a mapping matrix of weighting functions applied to 
the assumed basis functions. 

From previous arguments, we know that if a solution 
femoduation (54) exists and is unique for all gl. then 


Equation (54) can be transposed to: 
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In this form, the system of equations can be solved by 


imposing the boundary condition: 
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on ee, which then states that the tangential scattered field, 
Ee (which generates the current), is equal to the negative 


® s aL . e 
of the known impressed field, pe weteced tie cwlanren © onmn: 
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Hence, the system of equations can be solved for the exact 
current at the match points, and by enforcing current and 
charge continuity at segment ends, a very close approximation 
memene actual current distribution can be obtained. 

Equation (54) is solved by the Gauss-Doolittle method. 


Minis method the matrix [Z] is defined as two sub-matrices: 
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| = LJ (59) 


where sub-matrix LA] with elements Pas represents the electric 
field at the center of segment i due to the jth basis 


function, centered on segment j. Sub-matrix [BJ], with elements 


ae 


+40 represents magnetic fields generated by segment basis 


Lule 


Meetions. in the execution of this method, the matrix [Z] 
is factored into the product of an upper triangular matrix 


[U] and a lower triangular matrix [L]. That is: 


bat >= {L1 Lud Ss .2) 
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Equation (54) then becomes: 


Mei = CI) CLI) Cus (60) 


and LI] is computed in two stages: 


ma = (LJ CPI (61) 


end: 
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where Equation (61) is solved first for CF] (an intermediate 
Matrix) by forward substitution, and then for [I] by Lackward 
emesot il tution. 

Now that the current has been determined, the other 
non-zero scalar field components oy and Hs? can easily be 


computed. 
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Pole eeweoec er LON AND COMPUTATEZONAL RESULTS 


eeeeerODEL, DESCRIPTION 
The basic model used for this study is a simple monopole 
collinear with the Z-axis and mounted perpendicular to a per- 
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Fectly reflecting ground plane, that 


| 


fy plane. As can be seen from Figure 1, the antenna is 75 
meters in height and .6 meters in diameter. It is radiating 
at a frequency of 1 Megahertz ( the middle of the AM broad- 
cast band) with a wavelength (Az=c/f) of 300 meters; it is a 
quarter-wavelength monopole cperating at its resonant fre- 
quency. Transmit power is fed to the model threugh its base 
by a voltage that is normalized to obtain i000 watts of 
radiated power. 

Virtually all modeling schemes start from the basic 
premise that E- and H-fielid component values can be determined 
meom a knowledge of the current distribution on the antenna. 
Three generalized methods for describing tne current distri- 
bution will be discussed here: 1) assuming 4 current 
Miser ibuition, 2) modal modeling, and 3) sub-domain medeling. 

1. Assumed Current Distribution 

This is the classical method as described in [Ref. 4]. 
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the far-field, but in the near-field where stored energy 
peeongly depends on the current shape, large errors 
Gemerally occur. 
2. Modal Modeling 

To improve on the assumed distribution, the modal 
modeling method was devised. By this method the current 
distribution can be described as the sum of an infinite 
system of equations which precisely describe the current 
Ge@eeribution in the limit. This is the approach used by 
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Pemazano (Ref. 7]. Using an infinite series of Bessel-Fourier 
Integrals, which can be truncated to allow for desired 
accuracy, a system of closed form equations can be generated 
which lend themselves to computer solution. This method 
enforces boundary conditions (l.e. E. = 0) over the surface 
Of the antenna with precision (in the limit), but there is a 
practical limitation that must be imposed on the size of the 
system of equations that precludes an exact solution. If the 
current is expressed by a sum of harmonics, the total number 
of harmonics which must be used in the computations to match 
Beundary conditions is related to the ratio of dipole height 
to gap size. 
3. Sub-domain Modeling 

To avoid much of the mathematical complexity of modal 
modeling, sub-domain modeling divides a given structure into 
a concatenated series of sub-structures, or elements. The 
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approximated by geometric distributions (sin, cos, and linear 
constants) on the elements. This is the method used by NEC 
[Ref. 2]. The NEC code uses the scalar form of the integral 
equations to determine the current in each of the elements, 
then combines the results from each substructure to deduce the 


performance of the total structure. 


B. MODEL CONFIGURATIONS 

in this study three configurations of the model were 
assumed and values for the tangential, radial, and peak (time 
and spatial) fields were calculated for each. The three 
configurations were for a structure modeled with 5, 15, and 
25 sub-domains. In each of the configurations, a null value 
For the tangential E-field was imposed at the center point of 
each sub-domain (this point is called the match point). Since 
a closed form solvable function has not been identified which 
describes the current on the antenna, while maintaining a 
tangential E-field equal to zero, it has become necessary to 
approximate the current distribution. This is done in NEC 


Giomom [li.D. 


D 


through the Method of Moments és discussed in §$ 
Using this method the current can be approximated to an 
arbitrary desree of accuracy by increasing N. Fortunately, 
reasonable accuracy is obtained in most situations with a 


relatively small number cf segments, typically N = 1O/i. 
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fre oP eCrAT TIONS/OBSERVATIONS/ANOMALIES 

Mirleeseudy Ustlazed three configurations to determine two 
relevant facts: a) whether the calculated value of the 
near-rield was dependent on the number of segments used, 
and b) whether the size of the segment conformed to the 
guidelines used in the far-field. 

A non-analytic overview of the plots generated by this 
Seudy indicated that the calculated value of the near-field 
was independent of the number of segments into which the 
structure was divided. This was determined by overlaying 
plots of five and fifteen segments on a twenty-five segment 
plot. An example of one of these plots is shown in Figure 4. 
The spurious field emitted from the base of the antenna has 
Significantly different values. These variations were the 
result of using different values of input voltage, when the 
configuration was changed, to maintain a normalized radiated 
power output of 1900 watts.* The values along the antenna 
and at its top are essentially the same, for respective 
MecteiOns On the antenna, for all three configurations. 

The NEC manual [Ref. 2] recommends segment lengths (4) 
Jess than A/10 for normal modeling, and for critical regions 
aA < 4/20 is recommended. For a wavelength cf 300 meters 


as used in this study, segment lengths of A < 300/10 = 30 


“Radiated power denvended on the magnitude of the E- 
each segment, which in turn depended on the segment 
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Pret Sz vel ter 3, 15. and 23 segment models. 
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Figure 4 





meters, or A < 300/20 = 15 meters are suggested. Lengths used 
in each configuration were: A, = 70/5 =a Meters , ee (ayy kale 
= 5 meters, and A. = io) Zoe 2 Meters. Again, froma non-= 
analytic overview of plots similar to Figure +, it was 
determined that reasonable additional resolution was obtained 
in going from a five segment to a fifteen segment configura- 
mime (i.e, A = 15 meters to A = 5 meters), but the increased 


resolution in going from the fifteen segment to the twenty- 


DeRiceecrerGo A = cl meters) 


Five segment configuration (i.e. A 
was not justified by the increased computer resources it 
consumed. From this observation it could be concluded that 
smaller segment lengths are required for near-field calcu- 
Myetons, and that the lengths are on the order of 1/5 the 
value specified for the far-fields (i.e. A<A/50 for normal 
Situations and A<i/100 for critical areas), Additional 
meses and models should be run to medify, or confirm this 
observation. 

Intensity of the near-field is independent of the 


segment size, but smaller segment size 18 required to get a 


ity 


Mmesolution or how that intensity varied as a function o 


position along the antenna. Based on these observations, 


en 


determination was made to restrict all further discussion 
to the fifteen segment configuration to avoid duplication 
Sie errort . 

iweeeemeenrse Of Ene Investigacion, pious of the E- and 


i-fields as a function of distance along the G-axis were run 





maemeeie Dagse Of the monopole to its heignt, and also from 
Bem oase —©O twice the height of the monopole. The plots to 
just the height of the monopole were run to obtain increased 
Beeowuction as a function of position along the monopole. 
Peam@amation of the single and double height plots, side by 
side, revealed no increased resolution or additionel infor- 
Meeeton Was Obtained from a single height plot that was not 
available on the double height plot, so only plots to twice 
the antenna height are discussed in this section. 

Additionally, closely spaced plots within a segment were 
run to determine if there was any variation as a function of 
intra-segment position. This type of plot was run at the 
center of the antenna for: a) the end of a segment, b) the 
center of a segment, and c@) a position half-way between the 
positions deseribed ina)andb). This type of plot was also 
[ermeeear: a) the top of the antenna, b) a quarter of a segment 
above and below the tov, and ec) a haif a segment above and 
memew che top. In both cases, there was no observable 
M@eeererence in the a), b), or ec) plots. Consequently, they 
Maeve not been included in the discussion of plots that 
mollLows. 

Momameesule ct the above culling, the following are the 


plots included in Appendix A: 
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mee the height of the moncpole, at di 


ime 100), and 1000 ~adai fromethe axis. 
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me £-Field plots (Eos E ceees ne eLonsonr REO 


rae E peak 
a distance of one (1) wavelength from the monopole, for 
frees Of 2/100cths, 1/4, 1/2, 1, and 2 times the height 
of the monopole. 
feeeei-tield plots (H) as a function of Z to twice the 
mement Of the monopole, at distances of 1, 2, 5, 19, 100, 
and 1000 radii from the axis. 
feeeeei-rfield plots (H) as a funetion of R to a distance of 
See (l) wavelength for heights of 2/100ths, 1/4, 1/2, 1, 
and 2 times the height of the monopole. 
Fach of these categories are discussed in further detail 
below. 
ime ac Z) Plots 
he SiOWmeiner Leure | 5; E_fZ) along the antenna surface 
was much lower in magnitude (about 60 dB lower) than its value 
at the top. These lower values varied by an order of magni- 


ren 


tude as a result of NEC forcing E_ = 0 at match vboints and 


N 


approximating values (from the current distribution) else- 
where. These radical variations were quickly attenuated as 


the observation point (0.P.) moved away from the surface. At 


ais 
— Sed 


Hy 


ve radii (Figure §), the integrating effect 


H. 


meee stance 


O 


Ht 
—le 


}— 


of a wider field of view had essentially smoothed out the 
curve, and aiso decreased the difference in magnitude between 
values at the top of the antenna and those along its length. 
When the 0.P. moves to a distance of 100 radii, the integrating 


effect wilil neve created an almost uniform magnitude along 
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Plot Ez vs Z for a model of 15 segments. 


Plot Is at one radius. the Surface of the Monepole. 
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Plot Ez vs Z for a model of 15 segments. 


Plot Is of ten radil, 2.7 Meters above the Surface. 
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the antenna, and this uniform magnitude will continue to 
decrease in absolute vaiue as the distance of the 0.P. moves 
radially from the antenna. 

The magnitude of BE, and E. are approximately equal 
at the top of the monopole, but Ee maintains a value within 
an order of magnitude of its peak value all along the antenna 
fmmgmeciee (Figure 7). As a result of boundary conditions on £_ 
for the ground plane (similar to those for E on the antenna 
surface), the magnitude of EL quickly approacnes zero at the 
base. The shape of this curve remains pretty much the same 
as observations are made at greater distances. The magnitude 
of the curve decreases slowly and approaches a uniform value 
euonme the antenna. 
is a composite of Ee and E.: Near the antenna, 


la 
peak 


meeocs of E 
‘ peak 


implies that a is the dominant component in the near-field 


(Figure 8) strongly resemble plots of E This 


region. However, as the 0.P.'s are moved further outward, 
the E ‘E. resemblence fades and there is a peint around 
}/2r (see Figure 9) where the dominance begins to shift to 
E.. When E_ becomes completely dominant, the value of E ; 
Z, z peak 
is approaching a uniform value along the antenne. This 
results from the electromagnetic waves, generated by the 
point source at the top of the antenna, appreaching a plane 


Weave front at large distances. Alse nected in Figure 3 is the 


sharp decline of E_ above the height of tne antenna. This 


{\ 


indicates very little radiation in the axial direction, the 
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Plot Er vs Z for a mocel of 15 tagments. 


Plot 1s at one radius, the Surface of the Monopole. 
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Plot Epeok vs Z for a model of 15 segments. 
Plot Is ct one radius. the Surface of the Monopole. 
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Plot Epeak vs 2 fer a model of 15 segments. 


Plot ts ct 100 redii, 29.7 Meters above the Surtoce. 
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Peeacest percentage is directed radially outward. A more 
Pemprece set Of EC(Z) Plots are found in Figure A-1 through 
A-18 in Appendix A. 
mee a(R) Plots 

Figure 10 shows the very large value of ay Found at 
the top of the antenna decreases rapidly in magnitude as the 
O.P. moves radially away from the monopole. At a radial 
distance of approximately \/21the rapid decline in E, sub- 
sides and approximates the attenuation with distance 
observed at other points along the antenna (see Figure 11). 
Since Ee 0 at the surface, as shown in Figure ll, the 
Magnitude of the fieid rises rapidly with radial distance <o 
a peak at about 20 meters from the antenna, from which the 
magnitude slowly declines. 

mie ware OL decline of BS) appears to be uniform 
for all positions along the antenna, the only difference 


Deing the magnitude from which the decline begins. (Figure 12) 


Close to the antenna, the & mecusye GOoliows the £ 
peak r 
curve. At some unknown point, Ef) haceaim=nished in 
magnitude less than OR) 5 hence, the Eee curve switches to 
a 
a> low EfR) Gaieterer ts ).ca lie peingruae inten this switch 


Peemrs talls at a greater radiai distance for 0.P.'s at 


a larger starting 
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Sreater heights. This is the result o 


magnitude (at greater heights) from which E_ begins its 
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Unitorm rate of decline. For larger starting magntudes, a 
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Plot Ez vs r for co model of i5 segments. 


Plot is at the height of the monopscie along the r axis. 
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Plot Ez var for a model of 15 segments. 


Plot Js ot 1/2 the monopole heicht ciong the r axis. 
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Plot Er var for a modei of 15 tegments. 


Plot ts at 1/2 the monopole height along the r oxis. 
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Magnitude E,.,. (v/m) 


Plot Epeok ver for o model of 15 segments. 


Plot Is at the height of the monopole ciong the r cxis. 
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larger radiai distance is allowed for it to decline before 
the magnitude of E, exceeds a The point at which this 
occurs 1S so sharply defined, that the curve of E yeak has an 
feeparent discontinuity" in it. 


™7 


At twice the monopole height, the curve of Per 


(Figure 14) has an inverted cusp. This results from the same 
mrt in E/E, dominance discussed above. Directly above the 
antenna, Es. = 0; hence, a 1s larger than ae Puehas, he “Our. 
moves radially, E declines while EL rises to a peak. This 
Meem> the inverted cusp; then, EE, also begins a decline, but 


Meet aS rapid as Ee Mere plOtswors: 2 Jearemwound in Figures 
A-19 through A-33 in Appendix A. 
em GZ) Plots 


Pomecvescurrace of the monopole, the magnitude of H, 


is virtually constant; there is a slight decrease near the 

top (Figure 15). Above the monopole, there is a sharp drop 
moout SO dB's) in intensity of Hy Co ates LO tially 
Meee otene there. Similap plots of H, for greater distances 
from the antenna rollow the same type of curve, except, the 
magnitude of the field parallel to the antenna decreases with 
increasing distances from the antenna, and the sharp decline 
in intensity above the height of the antenna becomes less 


abrupt. At a distance of one wavelength, the difference 


between the magnitude of H. parallel to the antenna and above 


Te 


Meo height is about 20 dB. More plots of H(Z) are found in 
Pigure A-34 through A=-39 in Appendix A. 


iz 





Magnitude Eps ak (v/i) 


0.3 


Plot Epeak vs r for a mocel of 15 segments. 


Plot is at twice the monopole height along the r axis. 
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Plot Hoh ve Z to twice antenna height. 15 segments. 


Plot Is at one rodius, the Surface of the Monopole. 


Pirsure 15 
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Again, at the surface of the monopole the magnitude 
of His virtually constant. The magnitude of Eiccouccses 
rapidly with increasing radial distance, and uniformly for 
various antenna heights. Figure 16 is an example curve which 
delineates the benavior of H.all along the axis. As indicated 
above, for 0.P.'s directly above the antenna, Bue virtually 
non-existent. As the 0.P. moves radially outward, at twice 
the heignt of the antenna, HH pises quickly to a peak, then 
Slowly declines in magnitude. More plots of H(R) are 


found in Figure A-4G through A-44 in Appendix A. 
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Plot rphi vs r for a model of 15 segments. 


Piot is at 2/100th the monopole height clong the r axis. 
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ae eve See At TON OF NEC NEAR-FIELD CALCULATIONS 


Mee is still evolving. On-going efforts toward refining 
and improving its calculations are constantly occurring. In 
its present form, it has been extensively used, and docu- 
mented to provide accurate results in the far-field; its 
validity in close to the antenna is not so well documented. 
To show that NEC provides reasonable results in near-field 
regions, three tests have been run. They are described in 


the following paragraphs. 


eee OMPARTSON OF NEC WITH CLASSICAL ANALYSIS RESULTS 

Milis [Ref. 8], using equations similar to those found 
in Ref. 4 (pp. 333 to 338), calculated theoretical values 
OT L, as a function of radial distance from the antenna, 
for short antennas of various lengths. Two of his plots are 
Sie interest to this study; they are for divoles of half- 
Memeth (h) equal to .1A and .2\. The results of his 
calculations are plotted as broken lines in Figure 17, the 
plot of NEC's calculations for the same antennas are shown 
as solid lines. Remembering, the assumpticn of sinusoidal 


current is accurate only for the quarter-wavelength monovole®*, 


*As the antenna becomes shorter, the sinusoidal peak occurs 
Miiepqer back on the feed line and tne distribution on the 
antenna becomes more linear. 
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t 1s interesting to note NEC's calculations are more closely 


b 


matched as the antenna length approaches the resonant 
condition, especially for radial distances greater than 

A/2a (= 300/2 = 48 meters for this study). Balzano [Ref. 7] 
points out, that inside this distance, the classical apprcach 


".,.because nowhere in the near field (p<i/27) 


breaks down, 
does the dipole look like an elementary source (infinitesimal 
dimension) from the observation point". Balzano also points - 
out that the classical approach has a constantly increasing 
a field up to the dipole surface which is incongruent with 
energy density boundary conditions at the antenna edge. 


Therefore, the leveling off of the magnitude of ES? as 


predicted by NEC, is the more reasonable plot. 


Bee COMPARISON OF NEC WITH EMPIRICAL ANALYSIS RESULTS 

Balzano in Ref. 7 (Part II) performed some measurements 
On a dipole constructed from bronze rods of .0015 meters 
Baametcer (= 3.75 X woe that was .188 meters in length 
Moe 712), including a .0018 meter gap at the center from 
mimech 1t was fed. 

A model of this antenna was constructed using NEC. The 
Bomet2Uraticon chosen was a monopole of 19 segments, each 
segment being .9049 meters 1lcng. The monopole was mounted 
above a perfectly conducting ground plane in a manner 


mimeler £O that used for this study. Plots of E_ and EE. were 


c 
Ie 


mm@meaena ftumetion of 4 for radial distances of .003, .00c, 


es 





Mies 2, .09, -34, .05, and .06 meters from the dipole axis. 
Biewresults of these plots are shown in Appendix B. With 

the exception of a scaling factor, the information contained 
in these plots is the same as that discussed in Section IV.C. 
Beemene study. To further compare the results of NEC with 

the measured values, the measured peaks of the square of the 
Emeerric field intensity, for both axial and radial polariza- 
tions, at the radial distances indicated above, were extracted 
from plots in Ref. 7. These values were plotted in Figures 

18 and 19 along with the square of the veak of the fields 
extracted from the plots of the model run by NEC. As can be 
seen, very close correlation between the model and the 
measurements were obtained. The largest discrepancy was at 
.003 meters distance where the calculated value was inordinately 
higher than the measured value, relative to other calculated/ 
measured pairs. Balzano's calculated value was also quite a 
hit higher, which would indicate the measured value was a Dit 
low. Credence for this argument is supported by Balzano's 
statement, "...at closer distances, capacitive effects 


pecame predominant". 


C. COMPARISON OF NEC WITH MAXWELL'S EQUATIONS 
Chang [Ref. 6] points out that "...on the surface of the 
monopole the field component Hy should be proportional to the 


M@eeome ai ctrs DUtton 1(Z2)...'. lhe NEC caleulated value cf 


H was compared to the current flowing on the model. The 
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mesults are shown in Table I. 
(eA Bie) an 


COMPARISON OF NEC COMPUTED H-FIELD AND CURRENT AT THE 
VONOP OLE SURFACE 
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Quite good agreement iS maintained all along the antenna. 
The largest variation occurs at the end of the antenna ana 
that variation is still less than 1%. 


Chang uses the difference between: 
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#45 application to the data from the study at sample points 


yields Table II. 
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The correlation between tne estimated value and the calculated 
value 18 quite good above approximately 15 meters in height 

on the monopole. The wide variation below that point can be 
attributed to insufficient significant figures in tne values 
Or Hs used to estimate the difference. Also, NEC was 


computing the field emitted from the base cf the monopole, 


Bu 








wnich is net considered in the estimate of Maxwell's 
equations. Considering these two factors, the correlation 
Memsincprisingly good. 

By 1s expected, by Maxwell's equations, to be zero on the 
Barace Of the monopole. In NEC, this condition is only 
enforced at the match points. The best estimate of how well 
NEC approximates Ei =0 can be seen from the plot of E64) at 
the surface of the monopole, Figure 20. The magnitude of 
E, that is radiated can be found at the top of the monopole, 
memexpected. This value, from Figure 20, is nominally 10009 
v/m; at other points along the antenna, the nominal value can 
be approximated as 1 v/m. This drop of three orders of 


magnitude, or 60 dB, is a reasonable approximation of zero 


by virtually any engineering standard. 
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Plot Ez vs Z for a modal of 15 segmants. 


Plot Is af one radius. the Surface of the Monopole. 
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Figure 2 
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Vibe SeCNGiUSTONS AND RECOMMENDATIONS 


From the results of Section V, it can be seen that NEC is 
an effective analysis tool for near-field calculations and 
can be used and trusted to give meaningful results. The 
tool cannot be applied thoughtlessly; care must be exercised 
to ensure the model accurately represents the actual pDhysical 
structure. Zoning considerations must be factored into the 
procedure to ensure that sufficiently small segments are used 
to get the desired resolution, but not so many as to be 
wasteful of resources. Electrical parameters must be con- 
Sidered to ensure they represent the operaticnal environment 
into which the structure will be immersed. When these items 
are adequately considered, NEC will provide the desired 


results in a timely and cost effective manner 


Areas of further study could include expanded investi- 


CL, 


gations into the maximum segment length for near-fiel 
Bemeulations, investigation of the characteristics of the 
field emitted from the base insulator region, and effects 

on the near-field of using multiple parallel wires with cross 
eonnections in the manner in which broadcast towers are 


eeeually constructed. 


37 





Plot Ez vs Z for a model of 15 segments. 


Plot Is at one radius. the Surface of the Monopole. 
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Figure A-1 


Figures A-1 through A-6 are plots of aemae GunctLOM. Of 2 
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Magnitude = (v/m) 
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Plot Ex vs Z for a model of 15 segments. 


Plo? Is at two radil, .3 Meters above the Surtace. 
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Figure A-2 
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Magnitude E. (v/m) 
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Plot Ex vs Z for a model of 15 segments. 


Plot Is af five radil. 1.2 Meters above the Surface. 
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Magnitude E. (v/m) 
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Plot Ez vs Z for a model of 15 sagmenis. 


Plot is at ten rodi!, 2.7 Meters above the Surioce. 
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gure A-4 
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Magaitudo E. (v/m) 


Plat Ez v3 I for a zodai of 15 segments, 


Plot ts at 100 resi, 79.7 Motars above the Suriace. 
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Figure A-95 
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Piet &2 ve Z for a model of 13 segments. 


Plat 's at 1000 radii, 299.7 Meters above the Surtzce 
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Plot Er vs Z for a model of 15 segments. 


Plot Is at one radius. the Surface of the Monopole. 
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Figure A-/ 
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meeeures A-/ through A-12 are plots of E, as a function of z 
memeverious distances from the Z-axis. ~ 
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Magnitude E, (v/m) 


Plot Er vs Z for a model of 15 soqments. 


Plot is af two rodil. .3 Meters above the Surface. 
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Magnitude E (v/m) 


Plot Er vs Z for a model of 15 segments. 


Plot Is at five radil. 1.2 Meters above the Surface. 
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Figure A-9 
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Plot Er vs Z for a mode! of 15 segmente. 


Piet is af ten radii. 2.7 Meters above the Surface. 
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Figure A-10 
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Piet & vs Z for « model of 15 tegaents. 


Ptot ts at 100 reaciil, 77.7 A\eters above the Surince. 
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Figure A-1l 
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Magnitude E (v/w) 
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Piot Er vs Z for a mecal of 15 segaente. 


Pet is et 1009 red, 299.7 Meters above the furtace 
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Plot Epeak vs Z for a model of 15 segments. 


Plot Is at one radius. the Surface of the Monopole. 
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Figure A-13 


megures A-13 through A-18 are plots 
Pees tor Various distances from the 
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Plot Epeak vs Z for a model of 15 segments. 


Plot Is at two radil. .3 Meters above the Surfacs. 
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Plot Epeak vs Z for a model of 15 segments. 


Plot Is at five radil, 1.2 Meters above the Surface. 
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Magnitude Cosa (v/m) 
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Plot Epeok vs Z for a model of 15 segments. 


Plot is at ten radil. 2.7 Meters above the Surface. 
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Ficure A-16 
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Meagaitude Exe, (v/m) 


Piet cpeck vs Z for a model of 15 segments. 


Piet is at 100 reel, 29.7 Meters above the Surface. 
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Figure A-17 
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Plet Speak vs Z for a model of 15 segments. 


Ptet is ot 1000 rad. 299.7 Moters abeve the Surfcce 
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Plot Ex vs r for a model of 15 segments. 


Plot is at 2/100th the monopole height along the r axis. 
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Figure A-19 


Figures A~-19 through A-23 are plots of aS Mewar Ener Oon gor hk 
for varicus heights above the ground plane. 
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Plot Ez vs r for a model of 15 segments. 


Plot is at 1/4 the monopole height along the r axis. 
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Figure A-20 
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Magnitude E. (v/m) 
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Plot Ex var for a model of 15 segments. 


Plot Is at 1/2 the monopole height along the r axis. 
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Figure A-21 
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Magnitude E, (v/m) 


Plot Ex vs r for a model of |5 segments. 


Plot is at the height of the monopole along the r axis. 
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Plot Ez vi r for a model of 15 segments. 


Plot Is at twice the monopole height along the r axis. 
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Plot Er vs r for a model of 15 segments. 
Plot Is at 2/100th the monopole height along the ¢ axis. 
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Figure A-24 


Figures A-24 through A-28 are plots of EL as a function of R 
for various heights above the ground plane. 
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Plot Er vs r for a model af |5 segments. 


Plot Is ot 1/4 the monopole height along the r axis. 
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Figure A-25 
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Plot Er vs r for a model of 15 segments. 
Plot Is at 1/2 the monopole height along the r axis. 
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Figure A-26 
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Plot Er var for a modal of 15 segments. 


Plot Is at the height of the monopole clong the r axis. 
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Plot Er vs r for a model of 15 segments. 


Plot Is at twice the monopole height along the r axis. 
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Plot Epeok vs r for a model of 15 segments. 


Plot is at 2/100th the monopole height along the r axis. 
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Plot Epeak vs r for a model of 15 segments. 


Plot Is at 1/4 the monopole height along the r axis. 
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Figure A-30 
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Plot Epeak vs r for a model of 15 segments. 


Plot is at 1/2 the monopole helght along the r axis. 
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Figure A-31 
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Magnitude E,-,, (v/m) 


Plot Epeak vs r for a model of 15 segments. 


Plot Is at the height of the monopole along the r axis. 
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Figure A-32 
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Plot Epeak vs r for a model of 15 segments. 


Plot is at twice the monopole height along the r axis. 
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Figure A-33 
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Plot Hphi vs Z to twice antenna height. 15 segments. 


Plot Is at one radius, the Surfoce of the Monopole. 
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Figure A-34 


Figures A-34 through A-~393 are plots of am Sscraduiunction of 4 
meme various distances from the Z-axis. 
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Plot Hphi vs Z to twice antenna height. 15 segments. 


Plot is at two radil, .3 Meters cbove the Surfaces. 
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Plot Hphi vs Z to twice cntenna height. 15 segments. 


Plot Is ot five radii, L2 Meters above the Surfaces. 


Figure A-36 
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Plot Hphi vs Z to twice antenna height. 15 segments. 


Plot is at ten rodi!. 2.7 Meters above the Surface. 
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Figure A-3/7 
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Plot Hpk vs Z te twice entenna height, 15 segments. 


Piet ts at $00 rodil, 29.7 Meters aaovsa the Surface. 
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Figure A-38 
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Plot Hphi ve r for a model of 15 segments. 


Pilot Is of 2/100th the monopole height along the r axis. 
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Figure A-40 


Figures A-40 through A-44 are plots of H, as a function of R 
for varicus heights above the ground plane. 
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Plot Hphi vs r for a model of 15 segments. 


Plot is at 1/4 the monopole height along the r axis. 
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Plot Hphi vs r for a model of 15 segments. 


Plot ts at 1/2 the monopole height along the r axis. 
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Plot Hphi ve r for a model of 15 segments. 


Plot Is at the height of the monopole along the r axis. 
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Figure A-43 





Plot Hphi va r for a model of 15 segments. 


Plot !s at twice the monopole height along the r axis. 
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APPENDIX 8B 


Plot Ex ve z for a model of 19 segments. 


Plot is .3 cm frem monepole axis cieng the z axis. 
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Figures B-l through B-8 are plots or LE, as a ime Om Con. 2 


for various distances from the Z-axis. 
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Piet Ex vs z for « model of if segments. 


Plot is .6cm from the monepole axis along the z axis. 
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Plot Ex vs x for a model of 19 segments. 


Plot ts | cm from the monopole axis along the x axis. 
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Figure B-3 





Piet Ex vs x for a model of IF segments. 


Plot is 2 cm from the monopole axis ciong the z axis. 
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Plot Ex vs x for a model of !9 segments. 


Plot is J cm froe the monopole axis aleng the x axis. 


Magnitude E. (v/w 
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Pilot Ex ve z for a model of I9 segments. 


Plot 1 4 cm from the monopole axis aiong the x axis. 
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Figure B-6 
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Plot Ez vs z for a model of 19 segments. 


Piot is 5 cm from the monopole axis ciong the z exit. 


0.06 0.02 9.04 0.06 6.06 0.10 0.'2 
Distance Z (na) 


0 ee ~ 


te + 


tn * 
te 


ay 





Magnitude E. (v/m) 


Plot Ex vs z for a model of 19 segments. 


Piot ts 6 cm from the monopole exis along the z axis. 
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Figure 8-8 
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Plet Er vs x for a model of 19 segments. 


Pilot 's .3 cm from the menopele axis along the x axis. 
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Figure B-9 


Figures B-9 through B-16 are plots of EL, as a Give C1 Om Ou wt 
for various distances from the Z-axis. ~ 
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Plot Er vs z for a model of 19 segments. 


Plot |s .6cm from the menopeole axis along the & axis. 
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Figure B-10 
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Plot Er vs x for a model of I? segments. 


Plot ts | cm from the monepele axis aleng the z axis. 
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Figure B-1l 
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Plot Er vs z for a medel of I? segments. 


Plot is 2 cm from the monopole axis along the z exit. 
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Figure B-12 
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Plot Er ve x for a model of 19 segments. 


Plot Is 3 em from the menopeie exis along the x exit. 
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Plot Er vs x for a model of !9 segments. 


Plot !s 4 cm frem the monopole axis along the x axis. 
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Figure B-14 
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Plet Er vs x for @ model of 19 segments. 


Plot ts 5 cm from the menepole exis along the x azis. 
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Figure B-lo 


L46 





Plot Er vs s for o model of 19 segments. 


Plot ts 6 om from the monopole axis aleng the x axis. 
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Figure B-16 
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